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We consider finite linear or cyclic crystalline structures with molecular cells having nar- 
row pre-thieshold nuclear resonance. We prove that if the real part of such a nuclear res- 
onance lies within the energy band (the convex hull of the energy levels) of the crystalline 
structure arising of a separated molecular level, then there exist molecular crystalline states 
that decay exponentially in time and the decay rate pjj'"^ of these states in the main order is 
described by the formula T^^ =4^ where a is the value of the residue of the molecular 

in) 

channel transfer function at the nuclear resonance point and T)^ is the nuclear resonance 
width. 



I. INTRODUCTION 

Molecules are usually treated as purely Coulombic systems, while the strong interaction 
between their nuclear constituents is assumed to play a negligible role. However any Coulombic 
molecular level lying above the lower threshold of the nuclear subsystem, is embedded in the 
continuous spectrum of the nuclear sub-Hamiltonian. The coupling between the molecular and 
nuclear channels, hence, turns this level into a resonance (see, e. g., Refs. [2, 13, 18, iJ] 
and references cited therein). Of course, due to the wide Coulombic barrier between the nuclei 
and the short-range character of the nuclear interaction, this coupling, and thus the width of the 
resonance, which determines the fusion probability of the nuclear constituents of the molecule, 
is in general extremely small. 

However, as pointed out in Jsl 0], the situation may be rather different if the nuclear sub- 
system of a molecule has a sufficiently narrow near-threshold resonance. Examples of such 
nuclear systems may be read off from the data presented in [ 10]. Among them are even custom- 
ary systems like p p ^^O and p ^^O 1 1, 20], i. e., the nuclear constituents of the water molecule 
H2O or the hydroxyl ion OH with O being the isotope ^^O. For LID and H2O the influence of 
near-threshold nuclear resonances on the molecular properties has been studied in [jsl 0, 0] by 
estimating the overlap integrals between the corresponding molecular and nuclear wave func- 
tions. The best known example of such phenomena is the muon catalyzed fusion of deuteron 
and triton in the dtji molecule, where the near-threshold nuclear resonance ^He(3/2+) plays a 
decisive role [9]. 

Being motivated by the above special cases, we deal in the present work like in iIsLlZIi with 
a rather general model Hamiltonian related to the ones considered by Friedrichs in 111 ill . This 
Hamiltonian consists of a nuclear part, a molecular part with eigenvalues embedded in the con- 
tinuous spectrum of the nuclear part, and a weak coupling term which turns these unperturbed 
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eigenvalues into molecular resonances. Since the model is explicitly solvable, the mechanism 
of formation of the resonances becomes clearly visible. 

The following property pointed out in [3, 4] appears, in particular, as a general feature: if 
the nuclear channel itself has a narrow resonance with a position close to the molecular energy, 
then the width (the imaginary part) of the resulting molecular resonance is found to be inversely 
proportional to the nuclear width. In other words, a large increase of the decay rate of the 
molecular state, i. e. of the fusion probability, is observed in this case. Such a coincidence of 
nuclear and molecular energies is, of course, a very rare phenomenon in nature. 

Influence of the narrow pre-threshold resonances on the properties of infinite crystalline 
molecular structures was studied in [7]. In the present work we concentrate on more realistic 
finite crystals. A goal of this work is to show that the decay rate of a molecular state with 
the energy close to a near-threshold resonance may be considerably enhanced when arranging 
molecular clusters within a finite crystalline structure. The reason is that in such a configuration 
the original discrete molecular energy turns into a set of energy levels. That is, even if the 
position of the nuclear resonance differs from the original molecular level, it can get within 
this set. This allows for a fine tuning by exciting the crystalline structure to energies as close 
as possible to the energy of the nuclear resonance. We show that the lattice states, which 
correspond to such an initial choice of their quasimomentum distribution, decay exponentially 
with a rate which is again inversely proportional to the width of the nuclear resonance. 



II. TWO-CHANNEL MOLECULAR RESONANCE MODEL 

In this section we recall our main reasoning fslS] regarding an influence of a near-threshold 
nuclear resonance on the width of a molecular resonance in the case of a single molecule. 



A. Description of the model Hamiltonian 

Let us consider a two-channel Hilbert space = J^i © consisting of a nuclear Hilbert 
space J^i (channel 1) and a one-dimensional molecular space = C (channel 2). The ele- 
ments of are represented as vectors u = (^^^^ where ui G J^i and U2 G (u2 is simply 

a complex number). The inner product {u,v),^ = (mi,vi) -\-u2v2 in is naturally defined via 
the inner products (mi,vi) in Mi and M2V2 in ^2- 

As a Hamiltonian in we consider the 2x2 operator matrix 

where h\ is the (self-adjoint) "nuclear Hamiltonian" in and A2 G R a trial molecular energy. 
A vector b G ^ provides the coupling between the channels. It should be mentioned that the 
Hamiltonian dU) resembles one of the well known Friedrichs models [11]. 

If there is no coupling between the channels, i. e. for b = 0, the spectrum of A consists of 
the spectrum of hi and the additional eigenvalue A2. We assume that the continuous spectrum 
ac(/zi) of the Hamiltonian hi is not empty and that the eigenvalue A2 is embedded in Oc{hi). It 
is also assumed that A2 is not a threshold point of Oc{hi), and that this spectrum is absolutely 
continuous in a sufficiently wide neighborhood of A2. 



Perturbation of finite-lattice spectral levels by nearby nuclear resonances 



3 



A nontrivial coupling [b 7^ 0) between the channels will, in general, shift the eigenvalue 
A2 into an unphysical sheet of the energy plane. The resulting perturbed energy appears as 
a resonance, i.e., as a pole of the analytic (or, more precisely, meromorphic) continuation of 
the resolvent r{z) = {A — z)^^ taken between suitable states (see, e.g., |19]). In the present 
work we assume that such a continuation through the absolutely continuous spectrum of h\ in 
some neighborhood of A2 is possible at least for the matrix element {r\{z)b,b) of the resolvent 
r\{z) = {hi —z)^^. Then, from the explicit representation for the resolvent r{z) iSIP], one can 
easily see that the operator- valued function P2{A— z)^^\^^ admits meromorphic continuation 
to the same neighborhood, too. 

The poles of r(z) on the physical sheet are either due to zeros of the transfer function (see 



M2(z)=A2-z-/3(z) 

or due to poles of the resolvent ri (z) (see |S0])- The latter correspond to the discrete spectrum 
of the operator hi which may determine part of the point spectrum of A. This is true, in par- 
ticular, for the multiple eigenvalues of /ii. In any case it is obvious that the perturbation of the 
eigenvalue A2 only corresponds to solutions of the equation M2{z) = 0, i. e., of 

z = A2-/3(z). (2) 

This equation has no roots z with Imz 7^ on the physical sheet. Therefore, being eigenvalues 
of the self-adjoint operator A, they have to be real. Thus, Eq. ^ may have solutions only on 
the real axis and in the unphysical sheet(s) of the Riemann surface of the resolvent ri (z) ■ 

We start with a brief discussion of the case where the nuclear channel Hamiltonian hi gener- 
ates no resonances close to A2 in a domain 2i of the unphysical sheet which ajoins the physical 
sheet from below the cut. This assumption implies that for a wide set of unit vectors b = b/\\b\\ 
the quadratic form /3(z) = \\b\\^{ri{z)b,b) can be analytically continued in ^. Moreover, under 
certain smallness conditions for ||Z7||, Eq. Q is uniquely solvable ifisll in ^ providing in the 
main order (see, e. g., iflilfl^ ^ 

Z2 = X2-{ri{X2 + iO)b,b)+o{\\bf). (3) 

llz^ii^o 

p(2) 

The real and imaginary parts of the resonance Z2 = given by 

E^^ = A2-Re(ri(A2 + iO)Z7,Z7)+o(||Z7f), 

rf = 2\m{ri{X2 + m,b) +o{\\b\\')- (4) 



B. Perturbation of the molecular resonance by a nearby nuclear resonance 

p(i) 

Our main interest concerns the opposite case of a nuclear resonance zi=E\^' — > 

0, with a real part E^^^ close to A2. For the sake of simplicity we assume the corresponding pole 
of ri (z) to be of first order. Let the element b E J^i be such that the function /3 (z) admits an 
analytic continuation into a domain ^ which contains both points A2 and zi- This domain. 
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moreover, is assumed to belong to the unphysical sheet which adjoins the physical sheet along 
the upper rim of the cut. In & the function /3 (z), thus, can be written as 

p{z) = ^+rHz) (5) 

zi—z 

with P'^^^{z) being a holomorphic function. For a fixed "structure function" b = b/\\b\\ we have 
\a\ — Ca\\b\\'^ with a constant Ca determined by the residue of ri{z) at z = zi. Note that this 
residue is usually expressed in terms of resonance (Gamow) functions (see for example [ 17]). 
In fact, we assume that the resonance corresponds to an "almost eigenstate" of hi. That is, 

in principle a limiting procedure T^^^ ^ is possible so that the resonance turns into a usual 
eigenvalue with an eigenvector i^i e J^i. More precisely, we assume 

^ = ^^^^+0(1) as 4'^^0 (6) 

with = ^ 0. This can be achieved, e.g., if the Hamiltonian hi itself has 

a matrix representation of the form © and the resonance zi is generated by a separated one- 

lens 
Let 



dimensional channel. In such a case we would have C^^ = 1 (for details see Ref. ||4|], Sec. II). 



Rea > and Ima ^ Rea (7) 

and, for zE ^, 

\Im [5''^z)\>cs^\\bf and \r^z)\ <Cs^\\bf ■ 

with constants > and > 0. Furthermore, the coupling between the channels in the 
Hamiltonian (HI) is assumed to be so weak that 

ir^(z)|<c^||^f «r(,i^ while \a\ = CaWbf <t: (r^^^y . (8) 

It can be expected that these conditions are fulfilled in specific molecular systems even under 
the supposition that the nuclear width pj^'^ itself is very small. 

After inserting ^ for /3 (z), Eq. Q turns into the "quadratic" equation 

(A2-z)(zi-z)-a + (zi-z)/3'-'^s(z)=0 

which can be "solved", i. e., can be rewritten in form of two equations 




A2+zi-/3>-«g(z) , Ifh-zi-rHz'^^ 



+ a. (9) 

By Banach's Fixed Point Theorem, each of the equations ® has only one solution in the domain 
^. In case of the sign "— " we denote the root of © by Zmch in case of the sign "+" by Zmoi- 
According to (3] and the roots Znuci and Zmoi of Q are essentially given by 

a a 
A2-zi-rg(zi) h^i ^'^^ 

Zmoi = A2-r^(A2 + iO) + ^ filon + ^^^^ 

'^2-2l-p^(A2+lO) A2-Z1 
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From the second condition dH) follows 



<ti rjj''*. Consequently, this term provides 



in Znuci a very small perturbation of the initial nuclear resonance zi- As compared to pj^ 
it represents also in Zmo\ a very weak perturbation of the molecular energy A2. However, as 
compared to the result (|3l), valid in case of a missing nearby nuclear resonance, it can be rather 

large. In particular, if the molecular energy A2 coincides with the real part e]^^ of the nuclear 

im) r^'"^ 

resonance zi, then Zj^oi = Er^ ~ '~2~ ^^^^ 

(m) , Ima Kea 
EP^X2-2^ and rP=4^. (12) 

^R 

The width of the molecular resonance Zmoi the presence of a nearby nuclear resonance zi, 

thus, turns out to be inversely proportional to the nuclear width . 

The second inequality dU), chosen as a condition for \\b\\ reflects the fact that the "usual" 

molecular width r^^"* is much smaller than the width of a usual nuclear resonance F^^"*, 

QFf «c^(f(,'))'. (13) 

This can practically always be assumed for concrete molecules. 

Under condition ® the value of Q = |a|/||Z?|p differs from zero, Q > C > 0, as y'^^^ — > 
0. Therefore, in the presence of a narrow (T^^^ -C Ca/c^) nuclear resonance close to A2 the 

(m) (2) 

molecular width F)^ ' is much larger than the molecular width F]^ ' observed in absence of such 
a resonance. In fact, this ratio is determined by the large quotient ^^y^. 

^ R 

m. MOLECULAR RESONANCES IN A FINITE CRYSTALLINE LATTICE 

Let us assume that the "molecules" described by the Hamiltonian ^ are arranged in form 
of an finite one-dimensional linear (chain) crystalline structure. To describe such a crystal we 
introduce the lattice Hilbert space 

^ = © ^(0 (14) 
7=1 

representing an orthogonal sum of the Hilbert spaces associated with the individual cells 

Here the subspaces M'^'^ = ^ and J^2^ = = C are exactly the same ones as in Sec. Inland, 
thus, = . The elements of the total Hilbert space 5f are represented by the sequences 

( u^)\ 

u = (m*^^-* , u^'^^ , . . . , u^"^) with components u^'^ — J-x where Mj G J^fi and Uj G ^ = C. 



"2 



The inner product in ^ is defined by (m,v)^= (m^'sv^''')^(,). The subspaces ^1 = © ' 

7=1 «=i 
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and ^2 = © ' ^^'•^ ^ = 5^1 © represent respectively the pure nuclear and pure molecular 

i=l 

channels. 

In the present section we will first deal with the Hamiltonian H acting in according to 

(Hu)'^^^ =Wu^'-^'^ +Au^'^ +Wu^'+^\ i = 2,...,n-l (16) 

where only the interaction between neighboring cells is taken into account and the interaction 
operator W is chosen in the simplest form 

:) 

with w being a positive number. Such a choice of the interaction corresponds to the natural as- 
sumption that the cells interact between each other via the molecular states, while the direct in- 
teraction between nuclear constituents belonging to different cells is negligible. We assume that 
the closed interval [A2 — 2w, A2 + 2w] is totally embedded in the continuous spectrum Oc{hi) of 
hi and, moreover, that no thresholds of Oc{hi) belong to this interval. For the sake of simplicity 
we also assume that the interval belongs to the domain !^ introduced in Sec. Inland that for any 
II e [A2-2w,A2-l-2w] 

Im(ro(M±iO)S,S) ^0. (18) 
Obviously, the Hamiltonian ([T6b is a self-adjoint operator on the domain Dom(//) = © D^'^ 

i=l 

with d(') = Dom(/zi) © C. The resolvent R{z) = {H - of H possesses a natural block 
structure, = {R{j,k;z)}, j-,k= l,2,...,n. The blocks i?(7,fc;z) satisfy the equations 

WR{i-l,k-z) + {A-z)RU,k;z)+WR{i+l,k-z) = 5jkl, (19) 
j,k= 1,2,. 

where 5jk stands for the Kronecker delta and / for the identity operator in the Hilbert space 
of cells. Hereafter we assume Imz 7^ so that the value of z automatically belongs to the 
resolvent set of the operator H. The blocks R(j, k; z) themselves possess a 2 x 2 matrix structure, 
R{j,k;z) = {Rmn{j,k;z)}, m,n = 1,2, corresponding to the decomposition Jif = M'l © J^. 
The set of the sequences fk,k= 1 , 2, . . . , n, with the elements 



V « + 1 

Pk = ^. (21) 
n + 1 

forms an orthonormal basis for the (n-dimensional) Hilbert space of "-element sequences of 
the form {x\,X2-, ■ ■ ■ G C, j = 1, 2, . . . ,n. The Fourier transform 



{Fu){pk) = a/ ^ sin(pu) (22) 
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in reduces Eq. (fT9l) to 

{A- z)R{pk,Pk'\z) + 2co^PkWR{pk,Pk'\z) = dkk'I k,k' = 1,2,...,r, 



(23) 



and the numbers R{pk:Pk'\z) represent the matrix elements of the resolvent R{z) in this repre- 
sentation. From (l23t it immediately follows that 



R{Pk,Pk'',z) = G{pk;z)d, 



'kk'' 



where 



/ ri{z)b{-,b)r,{z) 

n{z) + - 



G{p;z) 



M2{p;z) 
{■Mniz) 



\ 



M2{p;z) 



n{z)b \ 

M2{p;z) 
1 

M2ip;z) / 



(24) 



(25) 



Here, the scalar function M2(p;z) reads 

M2{p;z) = A2-z + 2wcosp-/3(z) 



(26) 



The numbers pk given by ^TH represent the quasimomenta of the finite crystalline structure 
under consideration. 

Consider now the time evolution of the system described by the Hamiltonian H starting from 
a pure molecular state (p = (pi® (p2, \\(Pm\\ ^ '^m, 1,2, with (pi = and \\(p\\ = \\(p2\\ = 1- 
The probability to find the system at a time ? > in the molecular channel is given by 



^mol(<P,0 = l|P2e"'^>f, 



(27) 



where P2 is the orthogonal projection in ^ on the pure molecular subspace 5^2- Obviously, one 
can represent the time evolution operator exp{—iHt) in terms of the resolvent 7? (z) = (// — z)^\ 



exp{—iHt) 



2n\ 



dze 



-izt. 



H-z) 



(28) 



where the integration is performed along a counterclockwise contour 7 in the physical sheet 
encircling the spectrum of the Hamiltonian H. 

According to Eqs. (l24ll and (l26b the operator P2{H — z)^^ \c^^ acts in quasi-momentum rep- 
resentation as the multiplication operator. 



1 



(piiPk)- 



P2R{z)cp {pj,) = 

J M2{pk,z) 

Here ^{pk) stands for the values of the Fourier transform (|22|l of the vector 

<P2 = (<P2^^<pf^■■■,<P2"^)• 
Hence 



(29) 



2;ri 



(P2{Pk)J{Pk,t) 



(30) 
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with 

Hn.,)=idz . "P^'^" . (31) 

{ x,(pt)-z-m 

Repeating almost literally the analysis of Section III in Jvl] one finds that the asymptotics of the 
term J{pk,t) as ? ^ reads as follows 

J{pk,t) = exp{-izmoi(pyt)?} 



X 



a 



1 - ^ -2 + O e\pk, X2{pk) + iO) 

+ exp{-iz„uci(;?)?} (32) 



X 



a 



+ 0{£\pk,Zl)) 



+ £iPk,t); 



where 



£{pk,z) = ^ -. (33) 

[A2(w)-zi-/3-g(z)]2 

The function £{p,t) = 0{\\b\\^) is always small, \£{pj) \ <^ 1. By (flOb and dTTT) for the positions 
of the resonance poles we then obtain 

ZnudPk) = Zl- . . , (34) 

A2 + 2wcos;?fc-zi 

ZmoliPk) = A2 + 2wCOSj!?A:+T— . (35) 

A2 + 2w cos Pk — zi 

The asymptotics (l32t implies 



Pmoi{(P: = E \JiPk: 01' MPk) P = £ exp{-rj;") (pk) t} \cp2ipk) \^ + £it) (36) 

k=l k=l 

where 



^^R\pk) = -21mzmoi{Pk) = -2Im^— — ^ . (37) 

A2 + 2wcos;?A.-Zi 

The background term £{t) in is small for any t > 0, £{t) = 0{\\b\\^) and \£{t) | < 1. 

Further, let us assume that the number n of cells in the lattice is large and the real part E^ '^ 

of the nuclear resonance zi belongs to the interval [X2 — 2w, A2 + 2w], that is l^^^"* — A2I < 2w. 
Then, one can always prepare an initial molecular state (p which decays via the nuclear channel 

Re a . . __ 

with a rate close to 4— r-y (cf. formula UTl ). Indeed, under the assumption this maximum 

^ R 

is given by 

max r^, Md) = 4— 7--. 

o<p<K 

The most appropriate is the monochromatic molecular state (p with the only nonzero component 
(piiPko) associated with the quasimomentum pk^^ closest to 



Pmax = arccos 



E^^^ — A2 
2w 
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In particular, if the values of (pjiPk) nonzero only for quasimomenta restricted by 



cos pk- 



2w 



r(i) 
4w 



with some small 5 >0, then the width F^™'* given by the relation (l37t varies in an interval lying 

, , 1 4Rea , 4Rea 
approximately between — --— and — 

^ 1 + 52 A^) A^) 

^ R ^ R 

In a similar way one also treats a one-dimensional cyclic crystalline structure. In this case 
the Hilbert space $f is the same as in (fT4l) but the operator H reads 



{Hu)^^ =Wu'^'-^^ +Au^^ +Wu^'^^\ i = 2,...,n-l (38) 

If the intercellur interaction is still given by ( flTl . the only difference in the analysis will be the 
use of another complete orthonormal set in the space Z^. Instead of the sequences (EDI) one now 
employs the ortonormal sequences 

A(j) =^exp(ip^7), 7 = l,2,...,n (39) 

with quasimomenta pk given by 

Ink 

Pk = , k=l,2,...,n. (40) 

n 

After the Fourier transform (E^ with --^exp(ip^.7) the matrix elements of the resolvent {H — 

z)^^ again acquire the form (l24b . (l25t with the transfer function M2{p\z) given by (l26t . Hence, 
one concludes with formulas like in (l30b - (l37b and then observes that if the number of cells is 
large enough it is possible to prepare pure molecular states that decay with the rate close to 

/I Reg 

^ R 

In the same way one can also consider the finite two- and three-dimensional crystalline 
structures arranged of the molecular cells described by the Hamiltonian ([l]). If the cell has a 
sharp near-threshold nuclear resonance with energy embedded into the convex hull of the arising 
crystalline molecular levels, one again will find an enhancement of the decay rate for particular 
molecular states. As in the case of the one-dimensional lattices these molecular states should 
decay with the rate close to 4^^. 

^R 
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Bo3MymeHHe ypoBHeii KoneHHoro KpHCxajiJia 

6jIH3KOJ]eiKaiI|HMH fl/iepHMMH pe30HaHCaMH 



A. K. MoTOBHJiOB, B. B. BenaeB 

JIadopamopun meopemuuecKou (pusuKu, OMMH 
141980 ffydHa Mockobckou o6ji., Poccum 

Mbi HccjieflycM KOHCHHtie nHHefiHtie h uHKUHHCCKHe KpHcxannHHecKHe cxpyKTypti, sjieMeHxpaHtie 
MOJieKyjiapHbie ancHKH Koxopbix o5jiaflaiox ysKHM npe^noporoBbiM pesoHancoM. Mbi /joKasbisaeM, 
Hxo ecjiH BemecxBCHHaa nacxt sxoro pesoHanca HaxoflHxca BHyxpH BbinyKjiOH o6ojiohkh (nojiocbi) 
ypoBHefi 3HeprHH KpHcxanjia, B03HHKaiomHX H3 HeKoxoporo BbifleneHHoro MOJieKyjiapnoro ypoBHH, to 
cymecxByiox MoneKyjiapHbie coxosHHa KpHcxanjia, cicopocxb 3KcnoHeHiiHaiibHoro pacna^a Koxopbix 
r^*"^ B cxapmeM nopaflKe onHCbiaexca ^opMynofi pjj'"^ = r^e a - BenHHHHa BbiHexa xpaHccJjep- 

(JjyHKiiHH MOJieKyjiapHoro Kanana b tohkc aflepnoro pesonaHca, a pjj"' - innpHHa axoro pesoHanca. 



